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ABSTRACT 

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to 
correlate with the physical properties of their hosts, including the mass of the dark- 
matter halo. Using these clues as a starting point many studies have produced models 
that can explain phenomena like the quasar luminosity function. The shortcoming of 
this approach is that working models are not unique, and as a result it is not always 
clear what input physics is being constrained. Here we take a different approach. We 
identify critical parameters that describe the evolution of SMBHs at high redshift, 
and constrain their parameter space based on observations of high redshift quasars 
from the Sloan Digital Sky Survey. We find that the luminosity function taken in 
isolation is somewhat limited in its ability to constrain SMBH evolution due to some 
strong degeneracies. This explains the presence in the literature of a range of equally 
successful models based on different physical hypotheses. Including the constraint of 
the local SMBH to halo mass ratio breaks some of the degeneracies, and our results 
suggest halo masses at z ~ 4.8 of 1O 12 - 5±O 3 M (with 90% confidence), with a SMBH 
to halo mass ratio that decreases with time (> 99%). We also find a quasar luminosity 
to halo mass ratio that increases with halo mass (> 99%). These features need to be 
incorporated in all successful models of SMBH evolution. On the other hand current 
observations do not permit any conclusions regarding the evolution of quasar lifetime, 
or the SMBH occupation fraction in dark matter halos. 
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1 INTRODUCTION 

The Sloan Digital Sky Survey (SDSS) has discovered lu- 
minous quasars at redshifts as high as z ~ 6.4, i.e., when 
the universe was only a billion years old. The super-massive 
black-holes (SMBH) powering these quasars have been es- 
timated to have masses of ~ 10 9 solar masses. There is 
evidence that SMBHs contained a larger fraction of the 
host bulge mass at high redshift (Merloni, Rudnick & Di 
Matteo 2004; Shields, Salviander, Bonnin 2006; Peng et 
al. 2006). There are, however, several other questions re- 
garding the galaxies that host these high redshift quasars, 
and the physics that govern their evolution that have re- 
mained largely unanswered. Traditionally, investigation of 
the evolution of high redshift quasars has revolved around 
construction of models for the luminosity function which 
can be compared with observation. In this way successful 
models can be labeled as possible candidates to explain the 
formation and evolution of the quasars. An equally valu- 



able, though a rarer approach is to rule out potential mod- 
els because they are unable provide an explanation for the 
observed population. This approach is important because 
while a successful model can only suggest a mechanism as 
one possibility, a model which does not reproduce an ob- 
served phenomena may be definitely ruled out. 

Many successful models have been produced to ex- 
plain the high redshift luminosity function (e.g. Haiman & 
Loeb 1998; Haehnelt, Natarajan & Rees 1998; Kauffmann & 
Haehnelt 2000; Volonteri, Haardt & Madau 2003; Wyithe & 
Loeb 2003; Di Matteo, Croft, Springel & Hernquist 2003). 
Some authors advocate the view that the simpler a model 
the better because simple models can more clearly eluci- 
date the important underlying physics. Others would ar- 
gue that more detailed numerical models that attempt to 
capture more of the complex non-linear processes involved 
offer a more realistic and reliable description. Models for 
the luminosity function, therefore, take a variety of forms 
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and complexity. These range from analytic models in which 
simple physical prescriptions are combined with the Press- 
Schechter (1974) mass function and merger rates; to more 
complex semi-analytic models with merger trees and param- 
eterised laws to describe physical processes; through to nu- 
merical simulations with quasar activity included through a 
semi-analytic prescription. The published models of all types 
are broadly successful in their reproduction of high redshift 
number counts. In general this is because if one adds the 
number of free parameters to the number of physical as- 
sumptions, then one gets a number of parameters that is 
comparable to or smaller than the effective number of con- 
straints (which in practice includes only the slope of the 
luminosity function and the evolution of slope with redshift, 
in addition to the density normalisation). 

In this paper we highlight the fact that many differ- 
ent models are able to reproduce the data on the high red- 
shift quasar luminosity function, and note that this implies 
some degeneracy among input physical parameters or mech- 
anisms. Indeed, one such degeneracy was pointed out in an 
early model for high redshift quasar evolution (Haehnelt, 
Natarajan & Rees 1998). Assuming that quasar number 
counts trace the number of dark-matter halo hosts, Haehnelt 
et al. (1998) found that the luminosity function at z > 2 
could be equally well described using a constant SMBH to 
halo mass ratio combined with a constant lifetime; or using 
a mass dependent SMBH to halo mass ratio combined with 
a lifetime that decreased towards high redshift. The latter 
possibility is expected from a feedback scenario for SMBH 
growth (Silk & Rees 1998) . The degeneracy described above 
would imply that one cannot obtain strong evidence for feed- 
back using only the high redshift luminosity function data. 
More generally, these types of degeneracies mean that it is 
difficult to know which physical features of the various mod- 
els are being constrained by the data. Indeed we would like 
to know whether it is possible to conclude anything from 
published models other than the fact that the population 
of high redshift quasars can be reconciled with the ACDM 
cosmology. 

One feature that all models have in common is that they 
are based on the density and evolution of the dark-matter 
halo population. In semi-analytic models this population is 
generally modeled according to the Press-Schechter (1974) 
mass function (with extensions) . Similar results are obtained 
in cases where cosmological N-body codes are used because 
the Press-Schechter mass function provides an accurate re- 
flection of simulations over the mass-range of interest. In 
common with previous work, we use the Press-Schechter 
mass function as a starting point from which the quasar pop- 
ulation is generated. However we do not attempt to model 
physical processes such as SMBH growth and feedback etc., 
as is done in traditional semi-analytic modeling. Instead we 
take the more general approach of applying arbitrary param- 
eterisations of processes like the evolution of quasar lifetime, 
and the evolution of the relation between SMBH mass, halo 
mass and redshift. While it is true that we assign free pa- 
rameters to describe processes that we do not understand, 
the procedure for finding sets of these parameters which sat- 
isfy the data in a statistical sense is fundamentally different 
from the usual procedure of semi-analytic modeling which 
aims to find a single successful model. We are interested in 
the full range of allowed parameters as well as which sets 



of parameters can be excluded, rather than in a single set 
of parameters that is able to describe the data. In this way 
we hope to provide guidance for future attempts to build 
physical models of SMBH evolution. 

The optical quasar luminosity function shows a peak in 
its evolution at z ~ 2 — 3. At higher redshifts the quasar pop- 
ulation grows with time, and it is natural to relate the rise 
of the quasar population to the rise of the dark-matter halo 
population (e.g. Haehnelt, Natarajan & Rees 1998; Haiman 
& Loeb 1998; Volonteri, Haardt & Madau 2003; Wyithe & 
Loeb 2003). However near a redshift of z ~ 2 — 3 the non- 
linear mass scale moves from galactic scale to cluster scale. It 
is thought that the rapid fall in the density of bright quasars 
below z ~ 2 is due to a combination of a dwindling supply 
of cold gas at late times (Kauffmann & Haehnelt 2000) with 
an injection of feedback from the quasars into the surround- 
ing IGM that prevents further gas accretion onto collapsing 
systems (Scannapiecco & Oh 2004). Following the peak of 
quasar evolution one can no-longer relate the growth of the 
quasar (or galaxy) population directly to that of the dark- 
matter halo population in any direct or model independent 
way. For these reasons we restrict our attention to quasars 
at redshifts beyond z ~ 3.7, where we can relate the evolu- 
tion in the quasar luminosity function directly to evolution 
of the Press-Schechter (1974) mass function. 

In § |5| we introduce our parameterisation of 
SMBH/quasar evolution, discuss the relation of these pa- 
rameters to quasar observables, and present our method of 
parameter estimation as well as our choice of prior proba- 
bilities for the different variables. Our approach is an ex- 
tension of the formalism described in Wyithe & Padman- 
abhan (2005, hereafter Paper-I). In JUand 21 we describe 
the constraints that the available data impose on physical 
and model parameters. The effects of scatter in the relation 
between quasar luminosity and halo mass and of the choice 
of parameterisation and prior probabilities are discussed in 
f|S]and Jjj] Finally we compare some simple physical models 
with our derived constraints (S0 before summarising our re- 
sults in SjS] Throughout the paper we adopt the set of cosmo- 
logical parameters determined by the Wilkinson Microwave 
Anisotropy Probe (WMAP, Spergel et al. 2003), namely mass 
density parameters of Q m = 0.27 in matter, Q,b = 0.044 in 
baryons, Qa = 0.73 in a cosmological constant, and a Hub- 
ble constant of Ho — 71 kms -1 Mpc -1 . For the primordial 
power-spectrum of density fluctuations, we adopt a power- 
law slope n=l, and the fitting formula to the exact transfer 
function of Bardeen et al. (1986). 



2 MODEL PARAMETERS AND 
CONSTRAINTS 

2.1 Parameters 

In this section we introduce our parameterisation of SMBH 
evolution. Our aim is to describe a parameterised, rather 
than a physical, model of SMBH evolution. We do not at- 
tempt to assign values for these parameters in order to rep- 
resent a particular physical theory. Rather our goal is to 
consider all values so as to elucidate the acceptable range of 
parameter values as well as any degeneracies. The param- 
eterisation chosen provides a general set of models within 
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which the various different possible physical models are con- 
tained. Therefore by using this set of models to constrain 
the allowed values for these parameters, we may determine 
which physical processes can be unambiguously extracted 
from the available data. 

We define halo mass M to be the mass of a halo hosting 
a quasar with a luminosity .M1450 = —26.7 at z = 4.8. 
The relation between quasar luminosity L and halo mass M 
may have a mass dependence and a redshift dependence in 
addition to a normalisation (Lq/Mo). We parameterise this 
dependence using two indices <5 and 7 as 



L(z 



(—)' 
L \M J 



(1 + *) 



(1) 



The Press-Schechter (1974) mass function n(M,z), 
[with the modification of Sheth & Tormen (2002) that will be 
adopted throughout our discussion] yields the number den- 
sity N(> M(z),z) of dark matter halos above some mass 
M(z) at redshift z. If luminous quasars reside in a fraction e 
of such dark-matter halos, then the observed number density 
of quasars is given by the product of two factors: 



N(< Ml4 50 ) = N(> M(z),z)t 
where 



(2) 



TEeminft/r 1 ^),!}, (3) 

t q is the (unknown) quasar lifetime and H~ 1 {z) is the Hub- 
ble time (see e.g., Estathiou & Rees 1998). Our fourth pa- 
rameter is introduced via the 2 dependence of t q . While e 
could also change with z due to various effects (dust obscu- 
ration, beaming angle etc.), we expect the dominant contri- 
bution to evolution in r to come from t q . Nevertheless we 
allow for growth of e with time. The evolution of r can be 
parameterised by 



r(£) 

TO 



(4) 



The parameters a, 7 and <5, in addition to log 10 M there- 
fore govern the evolution of SMBHs. In this paper our goal 
is to constrain the allowed values of these parameters using 
observations of high redshift quasars. (Note that the model 
described above assumes an arbitrarily chosen a power-law 
parameterisation. The effect of modifying this choice is in- 
vestigated in SjfJ) 

2.2 Observational constraints 

There are 3 observational constraints which we will use to 
identify the allowed range of parameters governing the evo- 
lution of SMBHs (a, 8 and 7) as well as halo mass M. Let 
us discuss key observations that constrain these. 

First, observations of high redshift quasars (Fan et 
al. 2001;2003;2004) reveal an exponential decline in the 
quasar population with redshift of the form 

*(Xi450 < -26.7, z) oc 10 s ° bs xz , (5) 

suggesting approximate constancy of B. As a measure of the 
rate at which luminous quasars appear, we therefore use the 
exponential slope (B) of tN(> M(z),z), which according 
to the above parameterisation is defined as 



B 



0.434- 



a , dlog 10 7V(> M,z) 



1 + z 



+ 



0.434 jd\og 10 N(> M,z) 



(l + z)S dlog 10 M 



(6) 



Second, observations of high redshift quasars (Fan et 
al. 2001) also reveal a power-law decline in the density of 
quasars per unit luminosity as a function of luminosity of 
the form 



(7) 



As a measure of the relative numbers of quasars with dif- 
ferent luminosities at a single redshift we use the power-law 
slope of n(M, z) 

1 ri log 10 (n) 



(3 



(8) 



<5dlog 10 (M)- 

The factor of 1/8 accounts for the nonlinear relation between 
luminosity and halo mass. 

Third, the ratio of SMBH to host dark-matter halo mass 
is determined locally. Given a SMBH to host dark-matter 
halo mass ratio R(M, z) at high redshift, we can compute an 
extrapolated ratio at z = for halos of mass Mo = 10 12 M© 
using the parameters 8 and 7. 

is 



Ro (lO 12 M ,O) =R(M,z) 



M 



1O 12 M 



(1 + 2)" 



(9) 



If we estimate the SMBH mass powering the high redshift 
quasars, then we can calculate R(M, z) for a given value of 
host mass M and compare the extrapolated value of Ro for 
halos of 10 12 A/q at z — with the observations. 

In this paper we assign Gaussian probabilities for the 
observed distributions of B, (5 and Ro- In paper-I, based on 
the data of Fan et al. (2001;2003;2004) we found that B ohs is 
well constrained with a mean B = —0.49 and variance AB — 
0.07. Fan et al. (2001) find /3 = -2.58 and A/3 = 0.23. Fer- 
rarese (2002) determines the ratio M bh /M for M = 1O 12 M 
and finds a significant dependence on the assumed mass- 
profile, with values ranging between log 10 Ro = —5 (for an 
NFW profile) and log 10 Ro = —5.6 (for an singular isother- 
mal profile). Since the true profile is likely to lie between 
these extremes, we adopt a mean of log 10 Ro = —5.3 and a 
variance Alog 10 Ro = 0.3. 

2.3 A-posteriori parameter estimation 

From Bayes theorem, the joint a-posteriori probability dis- 
tribution for the parameters of interest (a, 8, 7 and M) as 
well as as is 

d 5 P 



dad'ydSdlog 10 Mdag 



LbL^Lr 



dP ■ dP ■ dP 

prior prior p 



rior.obs 



dP n] 



d-, 



da 



dd 



dlog 10 M dag 



(10) 



In what follows we will consider both the prior and observed 
constraints for 8, hence the designation P pr ior,obs in equa- 
tion 11011 . Since as has an observed distribution, we may 
marginalise over its dependence 

d 4 P 



dad'ydSdlog 10 M 



dz 



das 



obs 
rfPpri 



d 5 P 



das dad'ydSdlog^ Mdas 



(11) 
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The likelihoods in equation 1101 for B and (3 are 

-1 , B-B S 2 ZLLt§-dl\2 

L B = e 2 t_ SB- J an d Lb = e 2 . 



(12) 



In what follows we will consider both the case in which 
the constraint of the SMBH-halo mass ratio is not included, 
(so that Lr = 1), and the case in which the constraint of the 
SMBH-halo mass ratio is included. In the later case (which 
requires an assumption about the SMBH mass (Mbh), and 
extrapolation of equation Q between z = and z ~ 4.8), 
we have 



L R = 



dlog 10 M b 



d log 10 M b 



-e 2 



H 1 



(13) 



where we have marginalised over a prior dependence on 
SMBH mass powering the high redshift quasars. 

2.4 Prior parameter distributions 

In this section we specify the prior probability distributions 
(Pprior) for parameters under investigation. In cases where 
parameters are independently constrained we present our 
representation of the observed distribution (P b s )- 

The value of B is independent of the magnitude of the 
unknown quasar lifetime. However in addition to the halo 
mass, the value of B depends on the form of the redshift 
evolution of r. The z dependence of t q is handled by using 
the parameter a\ t . The expected range of ait may be set by 
the following considerations. First, if the quasar lifetime is 
determined by the mass e-fold timescale of the SMBHs, then 
t q is independent of redshift, r oc 1/H~ 1 (z) and an « 3/2. 
Second, if the quasar lifetime is determined by the dynamical 
timescale at z, then t q ~ H~ 1 (z) making r independent of 
redshift and an — 0. (This is also true if t q > H~ 1 (z).) The 
evolution of the parameter r also depends on the evolution 
of £ a (1 + z) a '- . We expect the value of e to grow with 
time, hence a e < 0. Thus the evolution of r is parameterised 
according to r oc (1 + z) a , where a = an + a e . We take the 
prior probability for a to be flat over a wide range of values 

dP ■ 

prior 



da 



1 

Qmax a K 

where a n 



-2, a n 



(14) 



In addition, the value of B depends on the evolution 
of the relation between luminosity and halo mass (7). In 
paper-I we have shown that 7 > under the assumption 
of constant e, by rejecting the null-hypothesis that 7 = 0. 
There are two components of the parameter 7. First, the 
SMBH to halo mass ratio can vary with redshift (a depen- 
dence described by a parameter 7 ra tio). Second the accretion 
rate could vary with redshift (a dependence which we de- 
scribe by a parameter 7 acc ). Since the accretion rate is close 
to Eddington at high redshift (Fan et al. 2004), we expect 
7acc > 0. The combined evolution of the quasar luminosity 
to halo mass ratio may then be described by 7 = 7ratio+7acc- 
Feedback scenarios at constant accretion rate imply 7 = 2.5. 
We use the following broad prior-probability distribution for 
7 which includes these values. 



dP n _ 



where 



7n 



-l,7n 



9/2(15) 



d*y 7max 7min 

The value of j3 is independent of the unknown values 
of the quasar lifetime (including its evolution a) as well as 



7. However in addition to the halo mass, the value of /3 de- 
pends on the form of the relation between quasar luminosity 
and halo mass (8). As before there are two components to 

8. First, the SMBH to halo mass ratio can vary with halo 
mass at fixed redshift (a dependence described by a param- 
eter 5ratio). Second the accretion rate at fixed redshift could 
vary with halo mass (a dependence which we describe by 
a parameter S acc ). Since the accretion rate is close to Ed- 
dington in the most luminous high redshift quasars (Fan et 
al. 2004), we expect <5 a cc > 0. The combined evolution may 
then be described by 8 = <5 r atio + 5acc Locally, the SMBH 
to halo mass relation is steeper than linear (8 > 1), while 
feedback driven scenarios at constant accretion rate, as well 
as observation prefer 8 ~ 4/3 — 5/3. We consider two cases 
for prior knowledge of 8. First, we use the following prior, 
which ignores externally derived constraints 



dP ■ 
dS 



1 



where <5 min = 1/2, 5 max = 5/2.(16) 



Second, we assume that the accretion rate is insensitive to 
halo mass (<5 a cc = 0), so that at fixed redshift the SMBH 
mass depends only on the halo mass. Locally, direct esti- 
mates of SMBH and host mass can be made through ob- 
servations of galaxy dynamics. These observations reveal a 
correlation between SMBH mass and the characteristic ve- 
locity of the surrounding stellar spheroid (e.g. Merritt & 
Ferrarese 2001; Tremaine et al. 2002), and by extension of 
the host dark matter halo (Ferrarese 2003). We therefore 
also consider the case where 8 is hypothesised to be redshift 
independent so that the observed constraint can be applied 
at high redshift. In this case we take 



dP n] 



d8 



1 t 5-5 ^2 

■ 2 y-KX-> 



27rA<5 n 



(17) 



with a mean and variance of 8 = 1.5 and AS — 0.25. 

Comparison of the high redshift SMBH to halo mass ra- 
tio with the ratio expected by extrapolating the local value 
to high redshift requires an assumption about the SMBH 
mass. SMBHs powering the highest redshift quasars are be- 
lieved to have masses of Mbh ~ 10 9 AfQ (Fan et al. 2001). 
We consider the following prior probability for the SMBH 
mass (Mbh) powering the high redshift quasars. 

dP • 1 

l* J orior ^ 



dl0g 1Q Mbh log 10 Mbh, max — l0g 1Q M b h,min 

where log 10 M b h,min = 8.5, log 10 M b h,max = 9.5 (18) 

Of course, all constraints depend on the value of 
log 10 M. SMBHs powering the highest redshift quasars are 
though to have masses in excess of 10 s Mq, implying halo 
masses of M > 1O 9 M . Galaxies are not observed to have 
halos in excess of 1O 14 M . Hence we choose a prior proba- 
bility for log 10 M of the form 

nor J- 



dlog 10 M log 10 M max - log 10 Mmin ' 

where log 10 M min = 9, log 10 M max = 14. (19) 

Finally, The evaluation of the PS mass-function, and 
hence of the parameters B and j3 depend on the adopted 
cosmology, with the greatest dependence being on as- We 
utilise the observed distribution for as from WMAP (Spergel 
et al. 2003) 
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dP o] 



-(- 



(20) 



da 8 ^/27i : Acr 8 
with mean and variance of dg = 0.84 and Acts = 0.04. 

3 PARAMETER CONSTRAINTS 

In this section we present a-posteriori, joint probability dis- 
tributions for combinations of a, 7, S and log 10 M. These 
distributions were obtained by marginalising equation 
over the other parameters. For example the joint probability 
for 7 and log 10 M may be found from 



d 2 P 



djdlog 10 M 



da 



dS 



d 4 P 



dad'ydSdlog^ M 



(21) 



We also determine a-posteriori cumulative probability dis- 
tributions for individual parameters. For example, the cu- 
mulative probability for log 10 M is 



P{< log 10 M) oc 



da 



log 10 m 



dlog 10 M' 



d& I d-y- 



d 4 P 



dad-ydSdlog 10 M' 
We show results in 3 different cases: 



(22) 



• Case-I: We begin by determining the constraints that 
can be made considering only the limits on B and j3 derived 
from the SDSS quasars, i.e. Lr = 1 and assumption only of 
the prior probability (equation lKiH for S. In this case we are 
trying to determine 4 parameters from two constraints, and 
so are strongly under-constrained. 

• Case-II: We next determine the constraints that can be 
made considering the limits on Ro in addition to those on 
B and (3, thus Lr is given by equation 11311 . As in Case- 
I, we assume a flat prior probability (equation 1161 for 5. 
The addition of the third constraint partially breaks the 
degeneracies in Case-I. 

• Case-Ill: Finally we replace the prior probability (equa- 
tion EiJ f° r <5> with the distribution observed locally (equa- 
tion E2J- 



Results are presented in Figures 11131 for Cases- 
I— III respectively. The sets of panels in the up- 
two rows show contours of d 2 P / (dady)\ , 
1- 1 > : 1 J,, ,-'„.„;,,.,M oW d*P/(d 7 dlog 10 A/)|2 

For the 



per 
d 2 P/(dad5) 



obs' 



d 2 P/{d8d 1 ) 



d 2 P/(dddlog 10 M)\ ohs , and d 2 P/(dadlog 10 M)| ob3 
display of joint probability distributions, the solid, dashed, 
dot-dashed and dotted contours correspond to values that 
are 0.61, 0.26, 0.14 and 0.036 times the distributions peak 
value, whose position is marked by a dot. When applied to 
a Gaussian distribution, these values correspond to the 1,2, 
3 and 5-sigma levels respectively. We also show a-posteriori 
cumulative probability distributions for individual param- 
eters in the lower rows of Figures 11131 (dark lines) together 
with the corresponding prior probability distributions (light 
lines) . 

For each case we determine the values corresponding to 
the maximum likeli-hood solution normalised with respect 
to a maximum possible value of unity. Maximum likeli-hoods 
of order unity imply that the model is able to accommodate 
all constraints simultaneously. The values corresponding to 



the maximum likeli-hood solution are shown in Figures 1 1131 
(large crosses). For Case-I the maximum value of the like- 
lihood is given by Li = LqLb- For Case-II we include the 
constraint of local mass ratio so that the likeli-hood becomes 
Li 1 = L/sLbLr. Finally the addition of prior information on 
the parameter 8 leads to a maximum value for the likelihood 
in Case-Ill of L nl = L L B L Re - 62 /{2AS) . 
We consider each case in turn: 

Case-I: In the first case we determine what can be 
learned regarding the connection between SMBHs and host 
halos, using only the evolution in the number counts of 
quasars at high redshift and the slope of the high redshift 
quasar luminosity function. The results are shown in Fig- 
ure Q In this case we have two parameters that may be 
measured (/3 and B) but 4 model parameters that lack prior 
constraints. Clearly this problem is under-constrained. In- 
spection of equations JSJ and © shows that /3 depends only 
on M and 5, while B is a function of all four parameters. 
Thus we would expect a degeneracy in the determination of 
parameters S and M which is driven by the observed lumi- 
nosity function slope /3. This degeneracy is clearly seen in the 
center panel of Figure although it is broken at high and 
at low M by the additional constraints from B. Thus even 
in this under-constrained case we are able to obtain an esti- 
mate of mass, and show that <5 > 1. The preferred value for 
the parameter S is quite insensitive to the value of either a 
or 7. Because the measured value of B is degenerate between 
M (or S through the constraint of (3), 7 and a, there will be a 
degeneracy along a line in the three dimensional parameter 
space. As a result we see that no constraints are possible in 
the 7 — a plane, while there are strong degeneracies between 
a and M and between 7 and M. The maximum value of the 
likelihood is Li = 0.92, implying that the best fit model is 
a very good fit to the high redshift data. The set of best fit 
parameters is (log 10 M, 7, S, a) = (10.48,-1,0.84,-2). The 
best fit value of a is on the edge of the range of parameter 
space, which is indicative of the degeneracies between a and 
other parameters. 

Parameters for which the data are the restrictive have 
a-posteriori probability distributions that differ from the as- 
sumed prior-probability distributions. As a result the lack of 
restriction imposed on models for SMBH growth by the high 
redshift luminosity function may be clearly seen in the lower 
row of Figure While the high redshift luminosity function 
places limits on each of the parameters 5 and M, no limits 
may be placed on the parameters a or 7. Additional observ- 
ables must therefore be included in order to constrain the 
models. We next turn to inclusion of the likelihood based on 
the local SMBH to halo mass ratio. 

Case-II: We can break part of the degeneracy seen in 
Case-I by including the likeli-hood for the local SMBH to 
halo mass ratio (equation I13H as an additional constraint. 
This inclusion has a two-fold effect on the parameter fit. 
Firstly, it couples the dependence of the overall likeli-hood 
of a solution on 8 and 7 in a direct way. The likelihood for 
the ratio R will therefore be degenerate between S and 7. 
Secondly, the value of 7 is now constrained by the likeli- 
hood for R in addition to the likelihood for B. This breaks 
the three-way degeneracy between M, a and 7. The inter- 
play between the different likelihoods and parameters is now 
coupled so that simple interpretation of the solution is more 
difficult than in Case-I. However we now have 3 observables 
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Figure 1. Constraints on power-law parameters, Case-I: The various panels in the upper two rows show contours for the marginalised a- 
posteriori joint probability distributions for combinations of the parameters 7, 8, a and M. The dots show the location of the maximum 
likcli-hood within the marginalised distributions, while the cross shows the location of the maximum likeli-hood within the full 4- 
dimcnsional parameter space. The lower row shows marginalised a-posteriori cumulative probability distributions for 7, 8, a and M 
(dark lines). The light lines show the respective prior probability distributions. The large open circles in the upper panels, and the light 
dashed lines in the central and lower panels show the location in parameter space of models A and B described in S 171 



(B, j3 and R), from which to constrain 4 parameters. The 
solution must therefore still be subject to a certain degree 
of degeneracy. 

The results are shown in Figure 

H We find that the 
combinations of 7 with M, 8 with M and 7 with 8 are well 
constrained. In particular we find 7 > 0. This may also be 
seen by inspecting the a-posteriori cumulative probability 
distributions. We see that 7, 8 and log 10 M have a-posteriori 
distributions that differ substantially from their prior prob- 
abilities, indicating that the data is able to constrain their 
values. The value of 8 lies between 1.1 and 1.6 (90%), and 
is larger than unity in agreement with the locally observed 
prior probability distribution for 8 that will be included in 
Case-Ill. The mass is constrained to be larger than 10 12 ' 2 Mq 
and smaller than 10 12 ' 8 M Q (90%), while 1< 7 < 2.7 (90%). 
Unless the accretion rate is changing rapidly with redshift 
near z ~ 4.5, the positive value of 7 implies a SMBH to 
halo mass ratio that decreases with time. The main degen- 
eracies due to the remaining degree of freedom are between 
a and M (lower right), between a and 8 (upper center) 
and between a and 7 (upper left), although the dependen- 
cies on a are weak. The maximum value of the likelihood 



is Ln — 0.48, implying that the best fit model remains a 
very good fit to the data following the addition of the con- 
straint on the mass ratio. The set of best fit parameters is 
(log 10 M, 7, 8, a) = (12.3, 2.41, 1.19, -2.). The value of a re- 
mains at the edge of the range of parameter space, indicating 
that the degeneracies between a and other parameters is not 
removed through the inclusion of the mass ratio. 

Case-HI: Finally we include a prior probability for 8 
(equation 1171) based on observations of SMBHs at low red- 
shift. The results are plotted in the upper and central rows 
of figure |3] Strictly speaking the solution is still under- 
constrained since there are 4 parameters and three con- 
straints. However the observed prior distribution for 8 is 
sufficiently narrow that the hypothesis of the slope in the 
local relation (8) also holding at high redshift strongly con- 
strains the solution. The result is that both 7 and M are 
tightly constrained by the data. The maximum likeli-hood 
value is Lin ~ 0.40, with a set of best fit parameters 
(log 10 M, 7, 8, a) = (12.7,1.65,1.47,2.55). From the lower 
row of Figure |3] we see that only 7 and log 10 M have a- 
posteriori distributions that differ substantially from their 
prior probabilities. The data does not restrict 8 much be- 
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Figure 2. Constraints on power-law parameters, Case-II: The various panels in the upper two rows show contours for the marginalised a- 
posteriori joint probability distributions for combinations of the parameters 7, 8, a and M. The dots show the location of the maximum 
likcli-hood within the marginalised distributions, while the cross shows the location of the maximum likeli-hood within the full 4- 
dimcnsional parameter space. The lower row shows marginalised a-posteriori cumulative probability distributions for 7, 8, a and M 
(dark lines). The light lines show the respective prior probability distributions. The large open circles in the upper panels, and the light 
dashed lines in the central and lower panels show the location in parameter space of models A and B described in S 171 



yond the assumed prior probability in this case. However 
this indicates that the preferred value of 8 at z ~ 4.5 is 
similar to the local value as hypothesised. As in Case-II the 
data shows a definite preferred solution for these quantities. 
SMBHs form a smaller fraction of their hosts mass at early 
times than they do today, 1 < 7 < 2.8 (90%), with high red- 
shift host masses of 10 12 - 2_12 ' 8 M© (90%). The parameter a 
remains unconstrained. 



4 QUASAR LIFETIME 

If a fraction e of dark matter halos contain SMBHs, then 
the total lifetime of the quasar can be estimated (Martini & 
Weinberg 2001; Haiman & Hui 2001) by dividing the quasar 
density ^(Muso < — 26.7, z) by e times the density N(> 
M) of halos larger than M, and then multiplying by the 
Hubble time (for t q < H' 1 ), hence 



N(> M,z) 



(23) 



There is therefore a one to one correspondence between ei q 
and M. As a result the joint a-posteriori distribution, with 



M replaced by rf q may be obtained as 



d 4 P 

dad'ydSd log 10 (ei q ) 
HP ■ 

l*-* prio 



obs 



d 5 P 



J 8 das dad'ydSdlog-^Q 
where 



Mda s 



ob: 



rflogioM 
dlog 10 (rf q )' 1 > 



d l°gio( rf q) 



dlog 10 Mj qV ' 



(25) 



In the third row of figure [2] we show the joint proba- 
bility distributions (computed using equation 1241 for a and 
log 10 (e£ q ), 7 and log 10 (ef q ), and for 8 and log 10 (et q ). The 
distributions were computed using Case-Ill constraints. The 
behavior is similar to that for the distributions of log 10 M as 
expected. If the occupation fraction is of order unity, then 
we find lifetimes of 10 6-7 years, in agreement with estimates 
at lower redshifts. Observed quasar lifetimes therefore sug- 
gest an occupation fraction e that is of order unity. Since 
the Hubble time is ~ 10 9 years at z ~ 4.8, we find that the 
occupation fraction must be larger than e ~ 10 _3 -10 -2 . 
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Figure 3. Constraints on power-law parameters, Case-Ill: The upper and second rows show contours for the marginalised a-posteriori 
joint probability distributions for combinations of the parameters 7, 8, a and M. The third row shows distributions with M replaced by 
(etlt). The dots show the location of the maximum likeli-hood within the marginalised distributions, while the cross shows the location 
of the maximum likeli-hood within the full 4-dimcnsional parameter space. The lower row shows marginalised a-posteriori cumulative 
probability distributions for 7, S, a and M (dark lines). The light lines show the respective prior probability distributions. The large open 
circles in the upper panels, and the light dashed lines in the panels of the second, third and lower rows show the location in parameter 
space of models A and B described in S 171 
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Figure 4. Constraints on power-law parameters where scatter in the SMBH to halo mass ratio is included, Case-Ill: The various panels 
show contours for the marginalised a-posteriori joint probability distributions for combinations of the parameters 7, 6, a and M. The dots 
show the location of the maximum likeli-hood within the marginalised distributions, while the cross shows the location of the maximum 
likeli-hood within the full 4-dimensional parameter space. The large open circles in the upper panels, and the light dashed lines in the 
lower panels show the location in parameter space of models A and B described in §171 



5 SCATTER IN THE LUMINOSITY-HALO 
MASS RELATION 

Up till now we have assumed a one-to-one correspondence 
between halo mass and quasar luminosity (equation^. How- 
ever we would expect a scatter in this relation, due to (at 
least because of) scatter in the Mbh — M relation, and due 
to variation in the accretion rate. To assess the importance 
of this scatter, we make the following modification to the 
calculation already described. We have associated a quasar 
of luminosity L with a halo of mass M. We can regard the 
mass M as characteristic, and replace the density of halos 
n(M, z) with an effective density 

n aS {M,z) = —= / dm'n(M' ,z)e^- [m ^ )2 , (26) 

V27rAm J 

which is obtained as a weighted average of the densities 
of halos that house a quasar of luminosity L. Here m = 
log 10 M. Effective values of B, 7 etc. can then be calculated 
using n e ff(M, z). We have computed the resulting probabil- 
ity distributions for a, 7, S and log 10 M as before in or- 
der to assess the importance of scatter in equation (0. Dis- 
tributions corresponding to Case-Ill were computed taking 
Am = 0.5 (so that 68% of halos are contained within 1 
decade of mass), and are plotted in figure 2] 

The results are similar to those obtained in the ab- 
sence of scatter. In particular the contours for distribu- 
tions move by less than 1-sigma, and have quantitatively 
similar values. The maximum value of the likelihood is 
Lin — 0.38, implying that the best fit model remains a 
very good fit to the data following the addition of scatter in 



the SMBH-halo mass relation. The set of best fit parameters 
is (log 10 M, 7, 8, a) = (12.62, 1.84, 1.47, 3.59). 



6 DEPENDENCE OF RESULTS ON CHOICE 
OF PRIORS AND PARAMETERISATION 

In this paper our aim is to determine what physical proper- 
ties may be inferred from the high redshift quasar luminosity 
function within a general, rather than a specific framework. 
Two issues arise regarding the statistical robustness of our 
results. Firstly, having chosen a particular parameterisation 
for our model we then calculate likelihoods for parameters 
using that parameterisation (for example log 10 M, 7 and a 
in our prescription). We have chosen flat prior probability 
distributions for these parameters and hence find aposteri- 
ori probabilities for their values. Now if the observations are 
able to genuinely constrain the values of these parameters, 
then our results should be independent of the prior used 
(within reason). We would therefore like to know whether 
our results are sensitive to the choice of prior. Secondly, 
we have chosen a powerlaw parameterisation for the red- 
shift dependencies of the lifetime and SMBH to halo mass 
ratio. These parameterisations may or not provide an ade- 
quate description of reality. We would therefore like to know 
whether the choice of parameterisation affects our physical 
conclusions regarding the mass of host dark matter halos, or 
the evolution of the SMBH to halo mass ratio. We discuss 
these issues in the remainder of this section. 
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Figure 5. Constraints on power-law parameters assuming alternative priors, Case-II: The various panels show contours for the 
marginalised a-posteriori joint probability distributions for combinations of the parameters 7, <5, a and M. The dots show the loca- 
tion of the maximum likeli-hood within the marginalised distributions, while the cross shows the location of the maximum likeli-hood 
within the full 4-dimensional parameter space. 



6.1 Dependence on prior probabilities 

Up till now we have assumed prior probability distributions 
for dimensionless parameters (a, 7 and 8) that are flat (i.e. 
dfprior / da oc 1 etc). For the dimensioned quantity (M) we 
have assumed a prior probability that is flat in the loga- 
rithm [i.e. dP P rior/dlog 10 (M) oc 1]. While these are the nat- 
ural choices, we have repeated the calculations in Case-II 
assuming priors for 7 and 8 that are flat in the logarithm 
(for < 7 < 9/2 and 1/2 < 8 < 5/2), and a prior proba- 
bility dPpriov/dM oc 1, where 1O 9 M < M < 10 14 M Q for 
host halo mass. We have not changed the prior of the un- 
constrained parameter a. 

Figure |5] shows the resulting parameter constraints. 
The set of best fit parameters is (log 10 M, 7, 8, a) = 
(12.87, 1.08, 1.47,4.0). The constrained value for M is quite 
similar to the original case (figure |2Jl, thus the estimate of 
this (physical) parameter would appear to be insensitive to 
the choice of prior. The values of 8 and 7 also take simi- 
lar values indicating that they are constrained by the data 
rather than by the prior. The conclusion that the SMBH to 
halo mass ratio was larger in the past (7 > 0) is therefore 
not sensitive to the choice of prior. The value of a remains 
the least constrained, although positive values are preferred 
in this case. 



6.2 Exponential parameter constraints 

The parameters a and 7 describe redshift dependencies of 
the lifetime and SMBH to halo mass ratio respectively. Up 
till now we have specified power-law evolutions for these 
quantities, and constrained the allowed values of the cor- 
responding parameters. Of course it is only our theoretical 



prejudice that leads us to specify a power-law, since we have 
no observational evidence for this form of parameterisation. 
The assumption of a particular prior probability distribution 
for a parameter (Bayes Postulate) leads to prior probabilities 
for the physical parameter of interest (like the life-time or 
mass ratio) that depend on the parameterisation employed. 
We would therefore like to check that our choice of model 
does not qualitatively bias the results. To check the effect 
of our choice of parameterisation, we have re-calculated our 
results assuming exponential rather than power-law evolu- 
tion for those quantities that vary with redshift (i.e. r and 
L(z)/Lo). However we retain our power-law form for mass 
ratio at fixed redshift (8) since this is an observed rather 
than a postulated parameterisation. 

Thus equations Q and Q become 



L(z) _ Mbh(z) 



L 
and 



Mbh.o 



t(z) _ az 

TO 



(27) 



(28) 



which leads to evolution described by the equation 
91og 10 iV(> M,z) 



B = 0.434a + 
+ 0.4347 



dz 

dlog 10 N(> M,z) 



(29) 



dlog 10 M 

The constraint from the local mass ratio (equation |5J be- 
comes 

l-S 



Ro (10 12 M Q ,0) =R(M,z) 



M 



W 12 M Q 



(30) 
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Figure 6. Constraints on power-law (<5) and exponential (a and 7) parameters, Case-Ill: The various panels show contours for the 
marginalised a-posteriori joint probability distributions for combinations of the parameters 7, <5, a and M. The dots show the location 
of the maximum likeli-hood within the marginalised distributions, while the cross shows the location of the maximum likeli-hood within 
the full 4-dimcnsional parameter space. 



while the constraint based on the slope of the luminosity 
function (equation |HJ remains unchanged. Note that the re- 
lation of B to R as a function of 7 is altered by the new 
definition. 

Figure [(J shows the resulting parameter constraints as- 
suming the analysis described in Case-Ill. The constrained 
value for M is quite similar to the power-law case (figurc[3]l, 
thus the estimate of this (physical) parameter would ap- 
pear to be insensitive to the choice of parameterisation (as- 
suming a monotonic evolution over the range of interest). 
The maximum value of the likelihood is Lm — 0.43, im- 
plying that the best fit exponential model also provides a 
very good fit to the data. The set of best fit parameters is 
(log 10 M, 7, 8, a) = (12.7,0.66,1.4,-0.41). The values of a 
and 7 will clearly differ from the power-law case. However 
the conclusion that the SMBH to halo mass ratio was larger 
in the past (7 > 0) holds in the exponential case as in the 
power-law case. The exponential parameterisation allows a 
preferred value of a ~ 0, indicating an occupation fraction 
and/or duty-cycle that does not vary much with redshift. 



7 COMPARISON WITH SIMPLE PHYSICAL 
MODELS 

Before concluding, we apply our idea to a couple of sim- 
ple test case models. Among the successful models of the 
high redshift quasar luminosity function, two examples em- 
ploy particularly simple but different physical hypotheses. 
In the first, the quasar lifetime is constant (perhaps set 
by the Salpeter time) with an occupation fraction of unity, 
the SMBH makes up a constant fraction of the host halo 
mass, and the quasar shines at its Eddington rate (Haehnelt, 



Natarajan & Rees 1998; Haiman & Loeb 1998). We refer to 
this as model-A. In the second model, the SMBH mass is set 
by feedback (through energy conservation) over the dynam- 
ical time of the gas reservoir, the occupation fraction is of 
order unity and the quasar shines at its Eddington rate (Silk 
& Rees 1998; Haehnelt, Natarajan & Rees 1998; Wyithe & 
Loeb 2003). We refer to this as model-B. In the present con- 
text these two models may be described by the parameter 
sets (aA,7A,<5A) = (l-5,0,l) and (q b ,7b,5b)=(0,2.5,1.66) re- 
spectively. The locations of these models in parameter space 
are shown in Figures Ill5l by the the light dashed lines, and 
the letters A and B. 

Considering first the constraints in Case-I (Figure 0, 
we see that the degeneracies prohibit any conclusion regard- 
ing a or 7, while the preferred <5 lies between the two models. 
Indeed the two models may not be distinguished at all based 
on the luminosity function data alone. On the other hand 
the simple models described should also be able to accom- 
modate the constraint from the local SMBH to halo mass 
ratio. In Case-II (Figures and |SJ the different values of 
a for the two models are again indistinguishable. However 
in contrast the values of 7 and 8 are more consistent with 
the feedback driven model (model-B). Finally, if the prior 
probability for 8 is included (Case-Ill, Figures 13141 then we 
have ruled out model-A a-priori. However in this final case 
it is still instructive to investigate the consistency of model- 
B with observation. The model-B values of 8b = 1.66 and 
7b = 2.5 are consistent with the data. 

These results suggest that model-B is more consistent 
with the data than model-A. However it is interesting to 
note that the most likely values for the parameters 7 and 
8 lie between those of the two models, indicating that nei- 
ther describes the astrophysics in its totality. On the other 
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hand, Figure [3] demonstrates in a model independent way, 
that the data would prefer a SMBH to halo mass ratio that 
is larger at higher redshifts (as suggested by the evolution 
of the virial-mass — virial- velocity relation), combined with 
a SMBH to halo mass ratio that is larger in higher mass 
halos. These conclusions are insensitive to the addition of 
scatter in the luminosity-halo mass relation, to the form of 
parameterisation, and to the prior probabilities assumed. 
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8 SUMMARY 

In this paper we have investigated the constraints that ob- 
servations of the high redshift luminosity function place on 
SMBH evolution. Our approach differs from that tradition- 
ally employed. Rather than hypothesise a physical model to 
govern the formation and evolution of SMBHs, we have de- 
scribed a general parameterised model and determined the 
regions of the resulting parameter space that are admitted 
by the observations. The traditional approach allows one 
to demonstrate that a particular model is viable, but does 
not provide a test of its uniqueness, nor which aspects of 
the theory are constrained by the data. Conversely our ap- 
proach explicitly demonstrates the degeneracies among the 
evolution of different physical properties, and thus reveals 
which conclusions may be robustly drawn from the data. 

We find that in isolation, the available data for the lu- 
minosity function of high redshift quasars is unable to dis- 
tinguish between a range of different physical models due 
to a series of degeneracies between different physical prop- 
erties. The inclusion of the local SMBH to halo mass ra- 
tio as a constraint removes some of this degeneracy and al- 
lows us to determine the physical parameters that may most 
clearly be derived from the data. We find that the mass of 
the high redshift quasar host halos was M = 1O 12 ' 5±O ' 3 M 
(90%). The ratio of SMBH to halo mass was larger at high 
redshift with a dependence M bh oc M(l + 2: )( 1 - 9 ± - 9 ) [ or 
M hh oc Me (0 ' 7±0 ' 5)z ] (90%). In addition, the ratio of quasar 
luminosity to halo mass increases with halo mass. We find 
L oc M 1A±0 ' 25 (90%). In the instance of an accretion rate 
that is insensitive to halo mass this relation is in agree- 
ment with that observed in the local universe. In hindsight 
the latter result might have been expected since the Press- 
Schechter mass function declines exponentially towards high 
redshift, while quasar number counts are observed to de- 
cline only as a power-law. These conclusions are not sensi- 
tive to scatter in the relation of SMBH to halo mass, or to 
the choices of parameterisation or prior probabilities, and so 
must be implicit in any physical model developed to describe 
SMBH evolution. However the evolution of lifetime and oc- 
cupation fraction are not currently probed by the data. 

In the future, the analysis described in this paper will 
be enhanced by additional constraints. For example, a larger 
sample of quasars would allow the slope of the luminosity 
function f3 and the evolution of quasar density B to be deter- 
mined over a range of redshifts. These additional constraints 
should break some of the remaining degeneracies and pro- 
vide us with constraints on the evolution of quasar lifetime 
to go beside the current constraints available for the host 
halo mass. 
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